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' Predicting urban growth is important for practical reasons, and also for the 

CN ' challenge it presents to theoretical frameworks for cluster dynamics Re- 

^ ' cently, the model of diffusion limited aggregation (DLA) [^J5[] has been applied 

^ ' to describe urban growth [0, and results in tree-like dendritic structures which 

^ ' have a core or "central business district" (CBD). The DLA model predicts that 

■ there exists only one large fractal cluster that is almost perfectly screened from 

' incoming "development units" (people, capital, resources, etc), so that almost 

£■ 

I ' all the cluster growth occurs in the extreme peripheral tips. Here we propose 

^ ' and develop an alternative model to DLA that describes the morphology and the 
^ ■ 

• area distribution of systems of cities, as well as the scaling of the urban perime- 

^ ' ter of individual cities. Our results agree both qualitatively and quantitatively 



with actual urban data. The resulting growth morphology can be understood 
in terms of the effects of interactions among the constituent units forming a 
urban region, and can be modeled using the correlated percolation model in the 
presence of a gradient. 

In the model we now develop we take into account the following points: 
(i) Urban data on the population density p(r) of actual urban systems are known to 
conform to the relation 0, p(r) = e~^^, where r is the radial distance from the CBD 
situated at the core, and A is the density gradient. Therefore, in our model the development 
units are positioned with an occupancy probability p{r) = p{r) that behaves in the same 



fashion as is known experimentally. 

(ii) We also take into account the fact that in actual urban systems, the development 
units are not positioned at random. Rather, there exist correlations arising from the fact that 
when a development unit is located in a given place, the probability of adjacent development 
units increases naturally — i.e., each site is not independently occupied by a development 
unit, but is occupied with a probability that depends on the occupancy of the neighborhood. 

In order to quantify these ideas, we consider the correlated percolation model . In 

the limit where correlations are so small as to be negligible [jl2|, a site at position r is occupied 



if the occupancy variable u{r) is smaller than the occupation probability p(r); the variables 
u{r) are uncorrelated random numbers. To introduce correlation among the variables, we 
convolute the uncorrelated variables u{r) with a kernel G{r), and define a new set of random 
numbers r]{r)= J G{r-r')u{r')dr\ The numbers r]{r) have long-range power-law correlations 
that decay as r~° if, e.g., we choose G{r) = 1/(1 + r^)"/^, where r = |r|. This choice assures 
that the correlations extend to the entire lattice. The Fourier transform of G{r), needed in 
order to convolute the variables in the Fourier space, is G{q) = Kp{q), where Kp{q) is the 
modified Bessel function of order f3 = {a — 2)/A. The assumption of power-law interactions 
is motivated by the fact that the "decision" for a development unit to be placed in a given 
location decays gradually with the distance from an occupied neighborhood. The correlations 
have the effect of agglomerating the units around a urban area, and they also quantify the 
interactions between clusters that are not physically connected. The correlation exponent a 
is the only parameter to be determined by empirical observations. 

To discuss the morphology of a system of cities generated in the present model, we show 
in Fig. |l| our simulations of correlated urban systems for a fixed value of the density gradient 
A, and for different degree of correlations quantified by the exponents a. We note the effect 
of correlations among the clusters. The strongly correlated system of Fig. |I]a compares better 
with actual urban settlements (i.e. see Fig. ^) than the weakly correlated system of Fig. 

In the simulated systems the largest city is situated in the core, which is regarded as the 
attractive center of the city, and is surrounded by small towns. All towns are nearly compact 



near their centers and become less compact near their boundaries, in quahtative agreement 
with empirical data on actual large cities such as Berlin, Paris, London, etc. (see i.e. Refs. 

GS)- 

Next we test the proposed model quantitatively. The urban boundary of the largest city 
is defined to be the perimeter of the cluster connected to the CBD. Since p{r) decreases as 
one moves away from the core, the probability that the largest cluster remains connected 
decreases with r. Hence Vf, the mean distance of the perimeter from the center, is determined 
by the value of r for which p(r) equals the percolation threshold — i.e. = pc, so Vf = 

A-Mn(lM) 0,11. 

The urban boundary in the model has the scaling properties of the external perimeter 
of a correlated percolation cluster in the presence of a gradient, which we calculate to have 
fractal dimension Df. — 1.3 for the uncorrelated case, and Df, ~ 1.4 for strong correlations 
(a — ^ 0). These values are consistent with actual urban data, for which values of De between 
1.2 and 1.4 are measured |^. Near the frontier and on length scales smaller than the width 
of the frontier, the largest cluster has fractal dimension df ^ 1.9. In contrast to the value of 
De, we find that the value of df is not affected by the long-range correlations parametrized 
by a. 

So far, we have shown how the correlations between the occupancy probabilities can 
account for the irregular morphology of the towns in a urban system. As can be seen in 
Fig. the towns surrounding a large city like Berlin are characterized with a wide range 
of sizes. We are interested in the laws that quantify the distribution N{A) of the area of 
the towns in such a system. Here, A is the area occupied by a given town. To this end, we 
calculate the actual distribution of the areas of the urban settlements around Berlin and 
London. We first digitize the empirical data of Fig. 4.1 of Refs. (which are shown in Fig. 
^), and Fig. 10.8 of Ref. [0]. Then, we count the number of towns with area A, putting the 
result in logarithmicaly spaced bins, and averaging over the size of the bin. The results for 
the distributions N{A) calculated in this way are shown in Fig. where we see that for 
both Berlin and London, N{A) follows a power-law. 

3 



This new result of a power law area distribution of towns, N{A), can be understood in 
the context of our model. Insight into this distribution can be developed by first noting that 
the small "towns" surrounding the largest city are all situated at distances r from the CBD 
such that p{r) < pc or r > r j . Therefore, we find N{A), the cumulative area distribution of 
clusters of area A, to be 

N{A) = r n{A,p) dp ~ A-^^+^l'^f^\ (1) 
Jo 

Here, n{A,p) ~ A~'^ g{A/Ao) is defined to be the average number of clusters containing A 
sites for a given p at a fixed distance r, and r = 1 + 2/df. Here, Ao{r) ~ \p{r) — Pc\~'^^^ 
corresponds to the maximum typical area of a "town" situated at a distance r from the 
CBD, while g{A/Ao) is a scaling function that decays rapidly (exponentially) for A > A^. 
The exponent v = z/(a) is defined by ^(r) ~ \p{r) —Pc\^'^, where ^(r) is the connectedness 
length that represents the mean linear extension of a cluster at a distance r from the CBD. 

In our numerical simulations we find a drastic increase of with the increase of the 
long-range correlations {a — > 0) (Fig. ^) ||TT|. The exponent z^(a) affects the area distribution 
of the small clusters around the CBD (Fig. ^c), as specified by Eq.(|I]), and can be used to 
quantify the degree of interaction between the CBD and the small surrounding towns. For 
instance, for a strongly correlated system of cities characterized by small a, i^(a) is large so 
that the area Ao{r) and the linear extension C,{r) of the towns will be large even for towns 
situated away from the CBD. This effect is observed in the correlated systems of cities of 
Fig. |. 

In Fig. ^ we also plot the power-law for the area distribution predicted by Eq. (|l]) along 
with the real data for Berlin and London. We find that the slopes of the plots for both cities 
are consistent with the prediction (dashed line) for the case of highly correlated systems. 
These results quantify the qualitative agreement between the morphology of actual urban 
areas and the strongly correlated urban systems obtained in our simulations. 

Before concluding, we discuss a generalization of our static model to describe the dy- 
namics of urban growth. Following the pioneering work of Clark |Q, the density gradient 
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X{t) was calculated as a function of the time for many urban areas around the world [[16 



Empirical studies show that the population density profile of cities presents a remarkable 
pattern of decline with time, which is quantified by the decrease of A(t) with time (see Fig. 
^ and Table 4 in Ref. [^) . In the context of our model, this fiattening pattern can be 
explained as follows. The model of percolation in a gradient can be related to a dynamical 
model of units diffusing from a central seed or core ||T^ , p!5[| . In this dynamical system, the 
units are allowed to diffuse on a two-dimensional lattice by hopping to nearest-neighbor 
positions. The density of units at the core remains constant: whenever a unit diffuses away 
from the core, it is replaced by a new unit. A well-defined diffusion front, defined as the 
boundary of the cluster of units that is linked to the central core, evolves in time. The 
diffusion front corresponds to the urban boundary of the central city. The static properties 
of the diffusion front of this system were found to be the same as those predicted by the 
gradient percolation model [ll4| , p!5| . Moreover, the dynamical model can explain the decrease 
of A(t) with time observed empirically. As the diffusion front situated around Vf moves away 
from the core, the city grows and the density gradient decreases since X{t) oc 1/tj. Thus, the 
dynamics in the model are quantified by a decreasing A(t), as occurs in actual urban areas. 
These considerations are tested in Fig. ^6, which shows our dynamical urban simulations of 
a strongly interacting system of cities characterized by a correlation exponent a = 0.05 for 
three different values of A. Qualitative agreement is observed between the morphology of the 
cities and towns of the actual data of Fig. and the simulations of Fig. 
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FIGURES 

FIG. 1. Simulations of urban systems for different degree of correlations. Here, the urban areas 
arc red, and the external perimeter or urban boundary of the largest cluster connected to the CBD 
is light green. In all the figures, we fix the value of the density gradient to be A = 0.009. (a) and 
(6) Two different examples of interactive systems of cities for correlation exponents a = 0.6 and 
a = 1.4, respectively. The development units are positioned with a probability that decays expo- 
nentially with the distance from the core. The units are located not randomly as in percolation, 
but rather in a correlated fashion depending of the neighboring occupied areas. The correlations 
are parametrized by the exponent a. The strongly correlated case corresponds to small a (a ^ 0). 
When a > d, where d is the spatial dimension of the substrate lattice (d = 2 in our case) , we recover 
the uncorrelated case. Notice the tendency to more compact clusters as we increase the degree of 
correlations (a — > 0). (c) As a zeroth order approximation, one might imagine the morphology 
predicted in the extreme limit whereby development units are positioned at random, rather than in 
the correlated way of Figs, la and 16. The results for this crude approximation of a non-interactive 
(uncorrelated) system of cities clearly display a drastically different morphology than found from 
data on real cities (such as shown in Fig. 2a). The non-interactive limit looks unrealistic in com- 
parison with real cities, for the lack of interactions creates a urban area characterized by many 
small towns spread loosely around the core. 

FIG. 2. Qualitative comparison between the actual urban data and the proposed model, (o) 
Three steps of the growth with time of Berlin and surrounding towns (from Ref. [13]). Data are 
shown for the years 1875, 1920, and 1945 (from top to bottom). The flattening of the density 
gradient is evident and corresponds to the decentralization of the urban area as the city grows. (6) 
Dynamical urban simulations of the proposed model. We fix the value of the correlation exponent 
to be a = 0.05 (strongly correlated case), and choose the occupancy probability p{r) to correspond 
to the data of Fig. 2a. 
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FIG. 3. (a) Log-log plot of the area distribution N{A) of the actual towns around Berlin and 
London. In the case of Berlin, the data are shown for the years 1920 and 1945 (corresponding to 
the last two panels in Fig. 2a), while the data of London are for the year 1981. A power-law is 
observed for the area distributions of both urban systems. The dotted line shows the predictions 
of our model for the uncorrected case (slope= 2.45), while the dashed line gives results for the 
strongly correlated case (slope= 2.06). Note that the area distributions for both cities agree much 
better with the strongly correlated case. The images have been digitized using an Apple Scanner, 
of resolution 150 dots per inche, and in this figure A denotes the number of sites covered by 
the development units for a given town. (6) Correlation length exponent i^(a) as a function of 
the correlation exponent a. Notice the drastic variation of u with the increase of the correlations 
(a ^ 0). This, in turn, corresponds to more compact clusters as can be observed in Fig. 1. (c) 
Log-log plot of the area distribution A'^(^) calculated for the present model for different degrees of 
correlation. From top to bottom, a = 0.2, a = 0.8, a = 1.4, and uncorrelated case. The linear fits 
correspond to the predictions of Eq. (1) using the values of z/(a) from Fig. 36, and df = 1.9. 

FIG. 4. Semi-log plot of the density profile p{r) for the three different stages in the growth 
of Berlin shown in Fig. 2a, corresponding to the years 1875, 1920, and 1945 from bottom to top. 
Least square fits yield the results A ~ 0.030, A 2± 0.012, and A ~ 0.009, respectively, showing the 
decrease of A with time. We used these density profiles in the dynamical simulations of Fig. 2b. 
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